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ABSTRACT 


This publication contains the preliminary agenda and summaries for the Third Annual JPL 
Airborne Geoscience Workshop, held at the Jet Propulsion Laboratory, Pasadena, California, on 
June 1-5, 1992. This main workshop is divided into three smaller workshops as follows: 

• The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) workshop, on 
June 1 and 2. The summaries for this workshop appear in Volume 1. 

• The Thermal Infrared Multispectral Scanner (TIMS) workshop, on June 3. The 
summaries for this workshop appear in Volume 2. 

• The Airborne Synthetic Aperture Radar (AIRSAR) workshop, on June 4 and 5. 
The summaries for this workshop appear in Volume 3. 
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Algorithm in Both Emission and Absorption 

Roger N. Clark, Gregg A. Swayze, Christopher Koch, and Cathy Ager 
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AGENDA (CONTINUED) 


third Annual JPL Airborne Geoscience workshop: 

Airborne Visible/Infrared imaging spectrometer 

(AVIRIS) 

1:30 p.m. Multiple Dataset Water-Quality Analyses in the Vicinity of an Ocean Wastewater 
Plume 

Michael Hamilton, Curtiss O. Davis, W. Joseph Rhea, and 
Jeannette van den Bosch 

2:00 p.m. MAC Europe 9 1 : Evaluation of AVIRIS, GER Imaging Spectrometry Data for 

the Land Application Testsite Oberpfaffenhofen 

F. Lehmann, R. Richter, H. Rothfuss, K, Werner, P. Hausknecht, A. Muller, and 
P. Strobl 

2:30 p.m. Using Endmembers in AVIRIS Images To Estimate Changes in Vegetative 
Biomass 

Milton O. Smith, John B. Adams, Susan L. Ustin, and Dar A. Roberts 
3:00 p.m. Break 

3:30 p.m. Atmospheric Correction of AVIRIS Data in Ocean Waters 
Gregory Terrie and Robert Arnone 

4:00 p.m. The 1991 AVIRIS/POLDER Experiment in Camargue, France 

F. Baret, C. Leprieur, S. Jacquemoud, V. Carrdre, X.F. Gu, M. Steven, 

V. Vanderbilt, J.F. Hanocq, S. Ustin, G. Rondeaux, C. Daughtry, L. Biehl, 

R, Pettigrew, D. Modro, H. Horoyan, T. Sarto, C. Despontin, and 
H. Razafindraibe 

4:30 p.m. AVIRIS Sensor and Ground Data System: Status and Plans 
Thomas Chrien and Earl Hansen 

5:00 p.m. Wrap up. 

5:30 p.m. End of AVIRIS Workshop. 

5:45 p.m. Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 
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AGENDA 


Third Annual JPL Airborne Geoscience workshop: 

THERMAL INFRARED MULTISPECTRAL SCANNER 
(TIMS) 

June 3, 1992 
Von Karman Auditorium 
Jet Propulsion Laboratory 
Pasadena, California 91109 

WEDNESDAY, JUNE 3, 1992 

7: 15 a.m. Shuttle bus departs Pasadena Ritz-Carlton Hotel for JPL. 

7:30 a.m. Registration and continental breakfast at JPL. 

8:00 a.m. Welcome 

8:30 a.m. TIMS Performance Evaluation Summary 

Bruce Spiering, G. Meeks, J. Anderson, S. Jaggi, and S. Kuo 

A Quantitative Analysis of TIMS Data Obtained on the Leaijet 23 at Various 

Altitudes 

S. Jaggi 

Analysis of TIMS Performance Subjected to Simulated Wind Blast 
S. Jaggi and S. Kuo 

9:30 a.m. Sensitivity of Blackbody Reference Panels to Wind Blast 
Gordon Hoover 

10:00 a.m. Break 

10:30 a.m. Preliminary Analysis of TIMS Performance on the ER-2 
S.J. Hook, VJ. Realmuto, and R.E. Alley 

1 1:00 a.m. Comparison of Preliminary Results From Airborne ASTER Simulator (AAS) 
With TIMS Data 

Yoshiaki Kannari, Franklin Mills, Hiroshi Watanabe, Teruya Ezaka, 

Tatsuhiko Narita, and Sheng-Huei Chang 

1 1:30 a.m. Simulation of ASTER Data Using AVIRIS Images 
Michael Abrams 

12:00 noon Lunch 

1 :00 p.m. Application of Split Window Technique to TIMS Data 

Tsuneo Matsunaga, Shuichi Rokugawa, and Yoshinori Ishii 

1:30 p.m. Atmospheric Corrections for TIMS Estimated Emittance 
T A. Warner and D.W. Levandowski 

2:00 p.m. An Algorithm for the Estimation of Bounds on the Emissivity and Temperatures 
From Thermal Multispectral Airborne Remotely Sensed Data 
S. Jaggi, D. Quattrochi, and R. Baskin 
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Agenda (continued) 


Third Annual JPL Airborne Geoscience workshop: 

Thermal Infrared multispectral Scanner 
(TIMS) 

2:30 p.m. Multi-Resolution Processing for Fractal Analysis of Airborne Remotely Sensed 
Data 

S. Jaggi, D. Quattrochi, and N. Lam 
3:00 p.m. Break 

3:30 p.m. Preliminary Analysis of Thermal-Infrared Multispectral Scanner Data of the Iron 

Hill, Colorado Carbonatite-Alkalic Rock Complex 
Lawrence C. Rowan, Kenneth Watson, and Susanne H. Miller 

4:00 p.m. The Use of TIMS for Mapping Different Pahoehoe Surfaces: Mauna Iki, Kilauea 
Scott K. Rowland 

4:30 p.m. Ejecta Patterns of Meteor Crater, Arizona Derived From the Linear Un-Mixing of 
TIMS Data and Laboratory Thermal Emission Spectra 
Michael S. Ramsey and Philip R. Christensen 

5:00 p.m. The Use of TIMS Data To Estimate the SO 2 Concentrations of Volcanic Plumes: 

A Case Study at Mount Etna, Sicily 
Vincent J. Realmuto 

5:30 p.m. End of TIMS Workshop. 

5:45 p.m. Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 
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Agenda 


third Annual JPL airborne Geoscience Workshop: 

Airborne Synthetic Aperture Radar 
(AIRSAR) 

June 4 and 5, 1992 
Von Karman Auditorium 
Jet Propulsion Laboratory 
Pasadena, California 91109 

THURSDAY, JUNE 4, 1992 

7:15 a.m. Shuttle bus departs Pasadena Ritz-Carlton Hotel for JPL. 

7:30 a.m. Registration and continental breakfast at JPL. 

8:00 a.m. Welcome 

8:30 a.m. The NASA/JPL Three-Frequency AIRSAR System 

J. van Zyl, R. Carande, Y. Low, T. Miller, and K. Wheeler 

9:00 a.m. A Snow Wetness Retrieval Algorithm for SAR 
Jiancheng Shi and Jeff Dozier 

9:30 a.m. Comparison of JPL- AIRSAR and DLR E-SAR Images from the MAC Europe ’91 
Campaign Over Testsite Oberpfaffenhofen: Frequency and Polarization 
Dependent Backscatter Variations From Agricultural Fields 
C. Schmullius and J. Nithack 

10:00 a.m. Break 

10:30 a.m. Monitoring Environmental State of Alaskan Forests With AIRSAR 

Kyle C. McDonald, JoBea Way, Eric Rignot, Cindy Williams, Les Viereck, and 
Phylis Adams 

11:00 a.m. Comparison of Modeled Backscatter With SAR Data at P-Band 
Yong Wang, Frank W. Davis, and John M. Melack 

1 1 :30 a.m. SAR Backscatter From Coniferous Forest Gaps 
John L. Day and Frank W. Davis 

12:00 noon Lunch 

1 :00 p.m. Panel Discussion on Future Emphasis of the AIRSAR System 
J. van Zyl, Moderator 

1 :30 p.m. Retrieval of Pine Forest Biomass Using JPL AIRSAR Data 

A. Beaudoin, T. Le Toan, F. Zagolski, C.C. Hsu, H.C. Han, and JA. Kong 

2:00 p.m. Characterization of Wetland, Forest, and Agricultural Ecosystems in Belize With 

Airborne Radar (AIRSAR) 

Kevin O. Pope, Jose Maria Rey-Benayas, and Jack F. Paris 

2:30 p.m. Strategies for Detection of Floodplain Inundation With Multi-Frequency 
Polarimetric SAR 
Laura L. Hess and John M. Melack 
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AGENDA (CONTINUED) 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 

Airborne Synthetic Aperture radar 
(AIRSAR) 

3:00 p.m. Break 

3:30 p.m. Supervised Fully Polarimetric Classification of the Black Forest Test Site. From 
MAESTRO 1 to MAC Europe 
G. De Grandi, C. Lavalle, H. De Groof, and A. Sieber 

4:00 p.m. Relating Multifrequency Radar Backscattering to Forest Biomass: Modeling and 
AIRSAR Measurement 
Guoqing Sun and K. Jon Ranson 

4:30 p.m. Poster Previews 

5 :00 p.m. Poster Previews 

5:30 p.m. End of session. 

5:45 p.m. Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 

6:30 p.m. Reception and poster sessions at the Pasadena Ritz-Carlton Hotel. 

9:00 p.m. Close of reception and poster sessions. 
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AGENDA 


Third Annual JPL Airborne Geoscience Workshop: 
poster Session 

Thursday, June 4, 1992 
6:30 to 9:00 p.m. 

Pasadena Ritz-Carlton Hotel 

1. Processing of AIRSAR Polarimetric Data for Soil Moisture Estimation Over 
Mahantango Watershed Area 

K. S. Rao 

2. Evaluation of Polarimetric SAR Parameters for Soil Moisture Retrieval 
Jiancheng Shi, Jakob J. van Zyl, and Edwin T. Engman 

3. Interaction Types and Their Like-Polarization Phase-Angle Difference Signatures 
Jack F. Paris 

4. Application of Modified VICAR/IBIS GIS to Analysis of July 1991 Flevoland 
AIRSAR Data 

L. Norikane, B. Broek, and A. Freeman 

5. Radar Analysis and Visualization Environment (RAVEN): Software for 
Polarimetric Radar Analysis 

K.S. Kierein-Young, AB. Lefkoff, and FA. Kruse 

6. Measuring Ocean Coherence Time With Dual-Baseline Interferometry 
Richard E. Carande 

7. A Bibliography of Global Change, Airborne Science, 1985-1991 
Edwin J. Sheffner and James G. Lawless 

8. ATTIRE (Analytical Tools for Thermal Infrared Engineering) — A Thermal 
Sensor Simulation Package 

S. Jaggi 

9. Kilauea Data Set Complied for Distribution on Compact Disc 

Lori Glaze, George Karas, Sonia Chernobiejf, Elsa Abbott, and Earnie Paylor 

10. The JPL Spectral Library 0.4 to 2.5 Micrometers 
Simon J. Hook, Cindy I. Grove, and Earnest D. Paylor II 

11. Lossless Compression of AVIRIS Data: Comparison of Methods and Instrument 
Constraints 

R.E. Roger, J.F. Arnold, M.C. Cavenor, andJA. Richards 

12. Simulation of AVHRR-K Band Ratios With AVIRIS 
Melanie A. Wetzel and Ronald M. Welch 

13. Soil Conservation Applications With C-Band SAR 
B. Brisco, RJ. Brown, J. Naunheimer, and D. Bedard 
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AGENDA 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 

THURSDAY POSTER SESSION (CONTINUED) 

14. Comparison of Edges Detected at Different Polarisations in MAESTRO Data 
Ronald G. Caves, Peter J. Harley, and Shaun Quegan 

15. Identification of Erosion Hazards in a Mediterranean Environment 
M. Altherr, J. Hill, and W. Mehl 
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Agenda 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 

Airborne Synthetic Aperture radar 
(AIRSAR) 

June 4 and 5, 1992 
Von Karman Auditorium 
Jet Propulsion Laboratory 
Pasadena, California 91109 

Friday, June 5, 1992 

7:15 a.m. Shuttle bus departs Pasadena Ritz-Carlton Hotel for JPL. 

7:30 a.m. Registration and continental breakfast at JPL. 

8:00 a.m. Oceanic Features Detected by SAR in the Mediterranean Sea During the MAC 
Europe ’91 Campaign 
Werner Alpers 

8:30 a.m. SAR Observations in the Gulf of Mexico 
David Sheres 

9:00 a.m. Investigation of AIRSAR Signatures of the Gulf Stream 

G.R. Valenzuela, JS. Lee, D.L. Schuler, G.O. Marmorino, F. Askari, K. Hoppel, 
JA.C. Kaiser, and W.C. Keller 

9:30 a.m. Mapping of Sea Bottom Topography 

C.J. Calkoen, GJ. Wensink, and G.H.F.M. Hesselmans 

10:00 a.m. Break 

10:30 a.m. Sea Bottom Topography Imaging With SAR 

M.WA. van derKooij, GJ. Wensink, andJ. Vogelzang 

1 1 :00 a.m. AIRSAR Surveys of Upper-Ocean Fronts Off California and Hawaii 
P. F lament 

11:30 a.m. Preliminary Results of Polarization Signatures for Glacial Moraines in the Mono 
Basin, Eastern Sierra Nevada 

Richard R. Forster, Andrew N. Fox, and Bryan Isacks 

12:00 noon Lunch 

1:00 p.m. Detecting Surface Roughness Effects on the Atmospheric Boundary Layer Via 
AIRSAR Data: A Field Experiment in Death Valley, California 
Dan G. Blumberg and Ronald Greeley 

1:30 p.m. Extraction of Quantitative Surface Characteristics From AIRSAR Data for Death 
Valley, California 
K.S. Kierein-Young and FA. Kruse 

2:00 p.m. The TOPSAR Interferometric Radar Topographic Mapping Instrument 
Howard A. Zebker, Spren N. Madsen, Jan Martin, Giovanni Alberti, 

Sergio Vetrella, and Alessandro Cucci 
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AGENDA (CONTINUED) 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 

Airborne synthetic aperture radar 
(AIRSAR) 

2:30 p.m. Evaluation of the TOPS AR Performance by Using Passive and Active Calibrators 
G. Alberti, A. Moccia.S. Ponte, and S.Vetrella 

3:00 p.m. Break 

3:30 p.m. Fitting a Three-Component Scattering Model to Polarimetric SAR Data 
A. Freeman and S. Durden 

4 00 p.m. Application of Symmetry Properties to Polarimetric Remote Sensing With JPL 
AIRSAR Data 

S.V. Nghiem, S.H. Yueh, R. Kwok, and F.K. Li 

4:30 p.m. External Calibration of Polarimetric Radar Images Using Distributed Targets 
Simon H. Yueh, S.V. Nghiem, and R. Kwok 

5:00 p.m. Wrap up. 

5:30 p.m. End of AIRSAR Workshop. 

5:45 p.m. Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 
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N94-16596 

TIMS PERFORMANCE EVALUATION SUMMARY 


BRUCE SPIERING, 6. MEEKS, J. ANDERSON, S. JAGGI, 

NASA/STENNIS SPACE CENTER 
STENNIS SPACE CENTER, MS 39529 


1 . PURPOSE 

The purpose of this effort was to investigate 
the TIMS operation under various conditions to 
determine what effect, if any, there is on the TIMS 
performance. This effort was undertaken because of 
reports from investigators that the TIMS is not 
performing as expected when flown on the ARC C-130. 
This effort was initiated with a meeting at SSC 
with representatives from SSC, JPL and ARC. 


S. KUO 


2 . PROBLEM DEFINITION 


The problem with the TIMS performance on the 
C-130 has been characterized in different ways. 
The primary problem seems to be as follows: Given a 
blackbody temperature reading from the blackbody 
controller, different digital count values can be 
obtained from the blackbody surface as measured by 
the TIMS, under different conditions encountered on 
the C-130. These differences result in incorrect 
surface temperatures computed from the data. 
Surface temperatures computed from TIMS data 
collected on the C-130 can be as much as 5C 
different than the ground truth data. 

3 . EXPERIMENTAL PLAN 


A series of experiments were planned which 
would systematically illuminate possible conditions 
which have a measurable effect on the TIMS 
performance. The hypothesis was that ambient 
temperature and wind blast on the surface of the 
TIMS blackbodies was affecting the TIMS performance 
on the C-130. 


The first experiment was a flight test with 
the TIMS installed on the SSC Learjet. The 
hypothesis was that ambient temperature and wind 
blast would affect the TIMS performance. Five 

flights at different altitudes were made over a 
large body of water. The data was distributed to 
three different investigators for analysis. A 
detailed analysis of the data is found in "A 
Quantitative Analysis of TIMS Data Obtained on the 
SSC/Learjet 23 at Various Altitudes. ",S. Jaggi, and 
"TIMS Performance Evaluation Flight Over Ross P. 
Barnett Reservoir in the Pre-Dawn Hours of 1 
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February 1992" , Palluconi,et.al. 

The second experiment involved simulating wind 
blast on the surface of the TIMS blackbodies while 
the TIMS was operated in the lab. The hypothesis 
was that wind blast on the surface of the 
blackbodies, at various speeds and orientations 
would affect the TIMS performance. A detailed 
analysis of the data is found in "Analysis of TIMS 
Performance Subjected to Simulated Wind 
Blast" , Jaggi , Kuo . 

A third experiment is planned if necessary. 
This experiment will be to define a flight 
experiment on the C-130. This experiment will 
include measurement of environmental conditions at 
the TIMS entrance aperture. 

4. FLIGHT TEST RESULTS SUMMARY 

Beginning at 3AM CST, five flights were flown 
in sequence at altitudes of 2Km, 6Km, 12Km, 6Kro, 
and 2Km over the Ross Barnett Reservoir near 
Jackson, MS. Radio-sonde data and water surface 
temperatures were acquired. Surface temperatures 
were acquired with both mercury thermometers and 
with a hand held radiometer. These measurements 
differed by approximately 2C over all the sample 
points . 

The TIMS data used in the initial analysis 
were the high and low temperature outputs from the 
blackbody controller, the high and low digital 
counts from the surface of the blackbodies, and the 
temperature outputs from the two ambient 
temperature monitors on the TIMS. One of the 
ambient temperature monitors was located in the 
TIMS entrance aperture near the high blackbody and 
the other was attached to the outside of the 
spectrometer housing. 

The four parameters obtained from the 
blackbodies were used to compute the slope and 
offset of the System Transfer Equation (STE) . The 
ambient temperature measurements were compared to 
the data from each altitude. Surface temperatures 
were computed from the video data independently by 
each investigator. Each investigator chose to use 
different techniques to compute surface 
temperature . 

The blackbody parameters and the ambient 
temperatures have been compared for the different 
altitudes. The computed temperatures have been 
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2 A QUANTITATIVE ANALYSIS OF TIMS DATA OBTAINED ON THE 
0 O LEARJET 23 AT VARIOUS ALTITUDES 

r- 

s. Jaggi 


N 94-16597 


Lockheed Engineering & Sciences Company 
Stennis Space Center, Mississippi 39529 


Su mma ry 

A series of TIMS data acquisition flights were conducted on the NASA Learjet 23 at 
different altitudes over a test site. The objective was to monitor the performance of the 
TIMS (its estimation of the brightness temperatures of the ground scene) with increasing 
altitude. The results do not show any significant correlation between the brightness 
temperatures and the altitude. The analysis indicates that the estimation of the 
temperatures is a function of the accuracy of the atmospheric correction used for each 
altitude. 

Five flights were flown in succession at altitudes of 2 km (low), 6 km (mid), 12 km 
(high), and then back again at 6 km and 2 km. The area selected was the Ross Barnett 
reservoir near Jackson, Mississippi. The mission was flown during the predawn hours of 
Feb. 1, 1992. Radiosonde data was collected for that duration to profile the characteristics 
of the atmosphere. Ground truth temperatures using thermometers and radiometers were 
also obtained over an area of the reservoir. The results of two independent runs of the 
radiometer data averaged 7.03 ± .70°C for the first run and 7.31 ± .88°C for the second 
run. Brightness temperatures were obtained for each channel for all flights. These were 
calculated by assuming the emissivity of the target to be 1.0 and computing the 
corresponding blackbody temperature given the upwelling ground radiance for each 
channel. Of particular interest were the temperatures corresponding to the region for 
which the ground truth data was available. Those regions were identified in each of the 
flights. The statistics of the raw video data, the corresponding radiance incident at the 
sensor, the upwelling ground radiance after it has been corrected for the atmosphere and 
finally the brightness temperatures were computed for those regions of each flight. 

The results of the brightness temperature of the region from channel 5 of the TIMS - 
the band least affected by the column of water and ozone - do not reveal a significant 
variation with altitude. The temperature varied from a mean value of 7.26 ± .42 °C at 2 
km, to 8.22 ± .45 °C at 6 km and to 8.35 ± .38 °C at 12 km. Considering the system 
accuracy of the TIMS sensor and the cumulative effect of the errors that can be introduced 
at various stages of the process of converting the raw data to brightness temperature, a 
variation of less than a degree between the temperatures obtained is not considered 
significant. The variation does not seem to be a function of the altitude as there is no 
continuous trend of change in temperature values from the low to the high altitudes. 
Rather, there seems to be a jump from the low altitude temperature to the mid and high 
altitude values, which are relatively close to each other, in spite of the high altitude being 
twice as high as the mid altitude. 

In this paper, the methodology used to obtain the brightness temperatures is described. 
The results of each stage of data conversion are illustrated. Also included are statistics for 
all the channels at each altitude and graphs for results of analysis on channel 5. 
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compared to the ground truth data for the different 
altitudes. The temperature differential between 
the TIMS data and the ground truth data does not 
exceed 2C. The blackbody parameters, ambient 
temperature measurements and different altitudes 
show no significant correlation. Therefore the 
original hypothesis has been found to be false, 
i.e., when the TIMS is operated on the SSC 
Lear jet, there are no measurable effects on TIMS 
performance due to ambient temperature or wind 
blast on the surface of the blackbodies. 

5. TIMS LABORATORY WIND ANALYSIS RESULTS SUMMARY 

The TIMS was operated in the lab with various 
configurations of wind directed at the surface of 
the blackbodies. At the time of this writing a 
maximum wind velocity of 40mph was achieved. Below 
approximately 30mph, very little effect was 
observed in TIMS performance. However, with wind 
speeds of 40mph, the TIMS continued to perform as 
designed. The data collected did not support the 
hypothesis that wind blast on the surface of the 
blackbodies under lab conditions adversely affected 
the TIMS performance. 

6. TIMS PERFORMANCE ANALYSIS RESULTS SUMMARY 

The data collected during these experiments 
indicates that the TIMS is a robust instrument. It 
can be operated under a wide range of conditions 
and still perform within its design goals. 

At this time, it has not yet been determined 
if the experiments described above are sufficient 
to predict the performance of the TIMS on the ARC 
C-130 . The third experiment has not yet been 
performed. It is believed that given a carefully 
defined and executed data acquisition mission, the 
TIMS can and will produce useful data when flown on 
the C-130. 
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Table 1 Digitized Raw Video 


Table 2 Radiance at Sensor 
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ANALYSIS OF TIMS PERFORMANCE SUBJECTED TO SIMULATED 

WIND BLAST. 

S. Jaggi S. Kuo 

Lockheed Engineering & Sciences Company 
Stennis Space Center, MS 39529 


Summary 

This report describes the results of the performance of the TIMS when it is subjected 
to various wind conditions in the laboratory. Various wind conditions were simulated 
using a 24" fan or combinations of air jet streams blowing toward either or both of the 
blackbody surfaces. The fan was used to simulate a large volume of air flow at moderate 
speeds (up to 30 mph). The small diameter air jets were used to probe TIMS system 
response in reaction to localized wind perturbations. The maximum nozzle speed of the 
air jet was 60 mph. A range of wind directions and speeds were set up in the laboratory 
during the test. The majority of the wind tests were conducted under ambient conditions 
with the room temperature fluctuating no more than 2 °C. The temperature of the high 
speed air jet was determined to be within 1 °C of the room temperature. TIMS response 
was recorded on analog tape. Additional thermistor readouts of the blackbody 
temperatures and thermocouple readout of the ambient temperature were recorded manually 
to be compared with the housekeeping data recorded on the tape. Additional tests were 
conducted under conditions of elevated and cooled room temperatures. The room 
temperature was varied between 19.5 to 25.5 °C in these tests. 

i 

The calibration parameters needed for quantitative analysis of TIMS data were First 
plotted on a scanline-by-scanline basis. These parameters are the low and high blackbody 
temperature readings as recorded by the TIMS and their corresponding digitized count 
values. Using these values, the system transfer equations were calculated. This equation 
allows us to compute the flux for any video count by computing the slope and intercept 
of the straight line that relates the flux to tine digital count. The actual video of the target 
(the lab floor in this case) was then compared with a simulated target. This simulated 
target was assumed to be a blackbody at emissivity of °.95 and the temperature was 
assumed to be at ambient temperature as recorded by the TIMS for each scanline. Using 
the slope and the intercept the flux corresponding to this target was converted into digital 
counts. The counts were observed to have a strong correlation with the actual video as 
recorded by the TIMS. 

The attached graphs describe the performance of the TIMS when compressed air is 
blown at each one of the blackbodies at different speeds. The effect of blowing a fan and 
changing the room temperature is also being analyzed. Results indicate that the TIMS 
system responds to variation in wind speed in real time and maintains the capability to 
produce accurate temperatures on a scan line basis. 
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SENSITIVITY OF BLACKBODY REFERENCE PANELS 
TO WIND BLAST 


Gordon Hoover 
Jet Propulsion Laboratory 
4800 Oak Grove Drive, Pasadena, California 91109 


As part of the effort at Stennis and JPL to 
discover the root causes of TIMS' temperature 
calibration errors, we undertook a series of 
experiments to measure the sensitivity of a heated 
plate to cooling by wind blast. 

We set up a powerful blower which was capable of 
generating a jet of wind in excess of 200 miles per 
hour. In the jet we mounted an electrically heated 
copper plate, one quarter of an inch thick and six 
inches square. The electrical heaters were capable of 
delivering a total power of 800 watts. The power to 
the heaters was feedback controlled with reference to a 
thermistor mounted on the back of the copper plate. 

The plate was mounted about 2.5 feet from the blower 
nozzle, at about 45 degrees to the direction of the 
jet. The jet was wide enough to wash the whole surface 
and its temperature at the plate was about 28 C. 

In order to simulate temperature differences which 
approximated flight conditions, we ran the plate at 
about 65 C, for a delta T of 37 C, 

The plate was instrumented with thermocouples in 
an attempt to measure the strength of temperature 
gradients within the plate. We placed thermocouples in 
holes drilled into the edge of the plate, in holes 
drilled from face to face, and surface mounted with a 
clamp. We found it quite difficult to measure the 
plate temperature in a way that we felt was insensitive 
to errors caused by the wind. For instance, merely 
clamping the thermocouple junction bead onto the 
surface with a small fiberglass tab led to errors of 
over five degrees due to heat flowing into the 
wind-cooled thermocouple leads. 

We observed apparent lateral gradients as well as 
depth gradients, but in the most extreme cases, the 
temperature differences within the copper metal 
amounted to no more than 1.0 C, 

If one considers the area covered by one of the 
100 watt heater elements, one can easily calculate the 
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maximum gradient that can be sustained with the heater 
at maximum output power. The power density was 5.2 
watts per square centimeter. The handbook value for 
the conductivity of copper was 3.8 watts/cm*degree (C) . 
After some calculating, I arrive at the maximum delta T 
of 0.82 C. 

1) delta T = Power*thickness/ (area*conductivity) 

Since full output power was not required to 
maintain the 37 C temperature difference, the maximum 
delta T inside the plate must be less than this 
estimate, maybe even as low as 0.1 or 0.2 C. 

The next step was to measure the gradient across 
the paint layer. Since paint has a much lower 
conductivity than copper, maybe it could sustain a 
gradient large enough to account for the errors. Two 
techniques were attempted. 

Method 1 

Method 1 required the comparison of the radiance 
of the wind blown target to a reference target 
unaffected by the wind. For this second target we used 
a sheet metal horn immersed in a stirred hot water 
bath. The temperature of the water was adjusted to 
approximate the anticipated temperature of the interior 
of the plate by reference to thermocouples fixed to the 
immersed horn and imbedded in the plate. The blower 
and plate heater were then turned on and the plate was 
allowed to reach a steady state. 

At this point, brightness temperatures were 
measured for the target and the comparison horn with a 
Barnes PRT5 Precision Radiation Thermometer and with an 
Omega radiation thermometer. 

As confirmation of the accuracy of the comparison, 
measurements were made of the two targets after the 
blower had been turned off and the plate had reached 
its new steady state. Under these conditions, the 
gradients through the plate and paint are minimized. 

The difference in brightness temperature between the 
plate and the immersed horn under these conditions are 
entirely due to the actual plate-horn difference and to 
the possible non-unity of the plate emissivity. 

Method 2 

The paint used on the heated copper plate was 
presumed to account for the greater part of the 
temperature drop during the operation of the blower. 

It is relatively easy to calculate the expected effect 
from increasing the thickness of the paint and to argue 
backwards to derive a conductivity for the paint from a 
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comparison of the brightness temperature of two patches 
of different thickness. 

Thickened paint patches two centimeters square 
were applied with multiple coats to areas of the plate 
which had been nearly isothermal in prior observations. 
The thicknesses were measured with a depth micrometer. 

The plate and the blower were then turned on and 
allowed to reach a steady state. The brightness 
temperature across the plates was observed using an 
Inframetrics Model 520 Thermal Camera and images were 
recorded on video tape. The Inframetrics camera is 
equipped with the ability to display the temperature 
profile generated by a single scan line. The sweep 
across the unequal patches of paint showed the 
brightness temperature differences clearly and this 
information permitted the temperature drop across the 
paint layer. 

Once the conductivity was known, we were able to 
calculate the temperature drop to be expected across a 
single thin layer of paint . This value could then be 
compared with the temperature drop from method 1 . 

One thing we discovered in the course of our 
experiments was that the single thickness coat of paint 
that appeared to eye to be solid black, fell short of 
being a perfect blackbody. The emissivity was 
estimated by a a number of methods, as follows: 

Under conditions of no wind, with the plate heated 
to a steady state, the brightness temperature of the 
surface was compared to the temperature of the interior 
of the plate. The gradient across the paint was 
assumed to be zero, attributing the whole drop observed 
to the departure of the plate from blackness. The 
radiance observed under these conditions is the sum of 
the thermally emitted radiation and the reflected 
radiation from the surroundings where the experiment 
was being conducted. We took this latter temperature 
to be the ambient air temperature. 

2) Total Radiance - emissivity * BB (T plate) 

+ (1 - emissivity) * BB (T air) 

2a) emissivity = (total rad. - BB(T air) 


(BB (T plate) - BB (T air)) 

In these formulas, the simplified form for the 
reflected term which derives partially from 
Kirchhoff's Law was adopted in the absence of detailed 
knowledge of the angular dependence of the incident 
radiation and the reflectivity of the plate. 
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1 . INTRODUCTION 


The Japanese Advanced Spaceborne Thermal Emission 
and Reflection radiometer ( ASTER ) , being developed for a 
NASA EOS-A satelite, will have 3 VNIR, 6 SWIR, and 5 TIR 
(8-12um) bands. An Airborne ASTER Simulator ( AAS ) has been 
developed for Japan Resources Observation System 
Organization (JAROS) by the Geophysical Environmental 
Research( GER ) Corp. to research surface temperature and 
emission features in the MWIR/TIR, to simulate ASTER’s 
TIR bands, and to study further possibility of MWIR/TIR 
bands . 

ASTER Simulator has 1 VNIR, 3 MWIR (3-5um), and 20 
(currently 24) TIR bands. We collected data over 3 sites 
-- Cuprite, Nevada; Long Valley/Mono Lake, California; and 
Death Valley, California --with simultaneous ground truth 
measurements. We present preliminary data collected by AAS 
for Cuprite, Nevada and compare AAS data with TIMS data. 

2. AIRBORNE ASTER SIMULATOR DATA 


The AAS data were acquired during the afternoon of 
Dec. 5, 1991. Each image is 495 lines of 544 pixels(512 
data, 16 heated reference black body, 16 ambient reference 
blackbody). Each pixel is signed 15 bit data. Table 1 
lists the approximate center wavelength and full-width at 
half- maximum ( FWHM ) for each channel based on preliminary 
laboratory tests. Additional information on the AAS can be 
found in Watanabe, et . al .( 1991 ) . The raw AAS data was 
converted to radiance by using the following equation: 

Rad d - (DN„ - DN W )/(DN C - DN w )*(Rad c - Rad w ) + Rad w 

The radiances for each reference blackbody were 
calculated based on the temperatures recorded by the 
operator at the beginning of each image and the center 
wavelength of the nominal bandpasses. The temperatures 
of the black bodies were assumed constant for the entire 
image. DN W was determined, separately for each scan line, 
by averaging the DN values of the heated blackbody on that 
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scan line. DN< was determined in a similar manner for the 
ambient black body. After conversion to radiances, the 
radiances were converted to 8-bit integers and processed 
to create panchromatic, false color, and decorrelat ion 
stretch images. The normalization method was used to 
determine temperature and emissivity for each pixel and 
channel . 

Performance data for AAS were estimated from 
approximate spectral tests at NASA Stennis Space Center 
and laboratory test at GER . Performance was also estimated 
by the Cuprite data which has good outcrops of quartzite, 
opal i zed rocks and carbonate rocks. 

3. COMPARISON OF AAS DATA WITH TIMS DATA 

TIMS data for the western Cuptite were acquired on 
Sept. 1,1 990. This TIMS data covered the same area of AAS 
CPWA2 data. The performances of TIMS were well known. 

Comparison of decorrelation stretch images from TIMS 
data with those of AAS data shows that the TIMS image of 
band 1-3*4 is similar to a part of the AAS image of band 
2-5-6 or band 2*4*6 or band 2*5*7 for 4 selected sites 
(Figure 1) . 

The decorrelation stretch images from the AAS are 
not as clean as the TIMS decorrela tion images. This is 
caused by several factors. First, the TIMS data was 
collected over much higher ground temperatures. Second, 
the bandwidths of the AAS are narrower than TIMS for most 
channels. Third, the decorrelation stretch amplifies the 
noise level in an image. 

This comparison provides some spectral performances 
and band characteristics of the current AAS . 

4. CONCLUSION 


Estimated center wavelengths of the AAS in its Dec. 

1991 configuration are listed in table 1. The FWHM 
bandwidth for the Dec. 1991 AAS was 0.2-0. 6pm, varying 
with channel. The Dec. 1991 FWHM for most channels is 
narrower than TIMS. Accordingly, the S/N ratio of the Dec. 

1991 AAS field survey data is less than TIMS Sept. 1990 
field survey data. The Dec. 1991 AAS NEdT , estimated from 
field survey data, is about 0.3-1. OK in the TIR bands. 

Following approximate calibration measurements at 
SSC in April 1992, GER will adjust and modify the AAS to 
Improve its performance before the Sept. 1992 data 
col lec t ions . 
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Table 1. Estimated Characteristics Of AAS For Dec. 1991 Survey Flight. 
Channel Detector Estimated Center Estimated 


1 

System 

TIR1 

2 

TIR1 

3 

TIR1 

4 

TIR1 

5 

TIR1 

6 

TIR1 

7 

TIR1 

8 

TIR1 

9 

TIR1 

10 

TIR1 

11 

TIR1 

12 

TIR1 

13 

TIR2 

14 

TIR2 

15 

TIR2 

16 

TIR2 

17 

TIR2 

18 

TIR2 

19 

TIR2 

20 

TIR2 

21 

TIR2 

22 

TIR2 

23 

TIR2 

24 

TIR2 

25 

MWIR 

26 

MWIR 

27 

MWIR 

28 

VNIR 


Wavelength ( Dec . 

8.2 

8.6 

8.9 

9.3 

9.6 

9.8 

9.9 
10.3 
10.5 
10.7 


9.8 

9.9 
10.1 
10.3 

10.5 
10.7 
11.0 

11 . 5 


3.0-3. 6 
3. 6-4. 2 
4. 4-5.0 
0.75-0.85 


1) FWHM( Dec . 91 ) 


0.6 

0.5 

0.5 

0.3 

0.25 

0.2 

0.25 

0.2 

0.2 

0.2 


0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.4 

0.3 


TIMS 



DECORRELATION STRETCH IMAGE 
BAND 134 



Figure 1. Decorrelation Stretch Image of TIMS Band 134 and Those 
Images of AAS. (A) Band 256 (B) Band 246 (C) Band 257 
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1. INTRODUCTION 

Absorptions by the atmosphere in thermal infrared region are mainly due to 
water vapor, carbon dioxide, and ozone. As the content of water vapor in the atmosphere 
greatly changes according to weather conditions, it is important to know its amount 
between the sensor and the ground for atmospheric corrections of TIMS data (i.e. 
radiosonde). On the other hand, various atmospheric correction techniques have already 
been developed for sea surface temperature estimations from satellites. Among such 
techniques, Split Window technique, now widely used for A VHRR (Advanced Very High 
Resolution Radiometer), uses no radiosonde or any kind of supplementary data but a 
difference between observed brightness temperatures in two channels for estimating 
atmospheric effects. 

In this paper, applications of Split Window technique to TIMS data are discussed 
because availability of atmospheric profile data is not clear when ASTER operates. After 
these theoretical discussions, the technique is experimentally applied to TIMS data at 
three ground targets and results are compared with atmospherically corrected data using 
LOWTRAN7 with radiosonde data. 

2. SPLIT WINDOW TECHNIQUE 

The theory of Split Window technique is based on equations of radiative transfer 
at two different wavelengths(for a detailed theory, see McMillin el al. 1984 ). To 
eliminate atmospheric terms from these equations, several approximations and 
assumptions(surface emissivities in two channels are unity, atmospheric absorptions are 
small and so on) are introduced, and equation (1) is derived: 

T| — T 2 = ~ (T s — T| ) 

1 . . (1) 
where T j and T 2 are observed brightness temperatures in channel 1 and 2, T s is surface 

brightness temperature, and kj and k 2 are absorption coefficients in channel 1 and 2. 

Equation (1) shows that the difference between surface and observed temperatures(T s -Tj) is 

proportional to the difference between observed temperatures in different channels^ -T 2 ). 

A coefficient of (T s -Tj) is usually determined by linear regressions of computer- 
simulated data (Deschamps et al. 1980) or buoy and satellite data set (Strong el aL 1984). 

In this aspect, this technique is an empirical one. As a result, a constant value is added 
to (1) because of approximation errors and (1) becomes equation (2), 

T s = T 1 + a(T 1 -T 2 ) + b or T s = a,T 1 + a 2 T 2 + b (2) 

3. TIMS CHANNEL COMBINATIONS FOR SPLIT WINDOW 
TECHNIQUE 
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To determine suitable channel combinations of TIMS for Split Window 
technique, brightness temperatures observed in TIMS channels were calculated using 
LOWTRAN7 for six model atmospheres, three surface temperatures for each model 
atmospheres, two looking angles(nadir and 45°), and the airplane altitude 4km. Results are 
shown in Fig. la and Fig. lb. Linear relationships between temperatures are almost 
satisfied for TIMS channel 5 and 6 combination(5/6), and combination 1/3. Non- 
linearities for combination 1/5, 3/4, and 4/5 probably come from differences between 
model atmospheres and atmospheric ozone. These figures suggest combinations 1/3 and 
5/6 are best combinations for Split Window technique. 

4. ESTIMATION OF LAND/WATER SURFACE TEMPERATURE 
USING TIMS DATA 

TIMS data used in this study were acquired over the Jasper Ridge area(between 
San Francisco and San Jose, California), at 1 1:30 a.m., August 30, 1990. Surface 
radiation temperatures were simultaneously measured at three locations(Searsvil!e Lake, 
rubber running track, and asphalt parking lot) and a radiosonde was launched from near the 
lake. TIMS data were converted into radiance and brightness temperature using internal 
reference data(air blast effects on reference surface temperature were corrected according to 
Schmugge et ai 1990). Radiance values were atmospherically corrected using 
LOWTRAN7 with radiosonde data and convened into brightness temperatures at the 
ground level. 

Coefficients of equation (2) were determined by LOWTRAN7 simulations. In 
this case, five relative humidity profiles based on Mid Latitude Summer model were used 
in the simulations instead of six model atmospheres. Simulation results are shown in 
Fig. 2, and coefficients were determined as follows (subscripts indicate TIMS channel): 

T S =T 3 + 1 ,705(T 3 -T j )-0.94 and T s =T 5 +3.238(T 5 -T 6 )+0.03 (3) 

Results of temperature estimations are shown in Table!. Temperatures 
atmospherically corrected with Mid Latitude Summer model are higher than with 
radiosonde data because the model overestimates the amount of water vapor. Estimated 
temperatures using Split Window technique agree with temperatures atmospherically 
corrected with radiosonde data within 2°C for Searsville Lake and the running track. But 
for the parking lot, the discrepancy is relatively large probably because of spectral 
variations of the surface emissivity. 

5. CONCLUSION 

Split Window technique can be applied to TIMS data if the channel combination 
is suitable. But equation (3) may not be used for other cases because coefficients of 
equation (2) depend on the type of the atmosphere and the altitude of the airplane. 
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b) 



Fig. 1. Results of LOWTRAN7 simulations: T s is surface brightness temperature 
andT,,T 3 ,T 4 ,Tj andT 6 are remotely observed brightness temperatures in TIMS Ch. 


1, 3, 4, 5 and 6, respectively. 



Fig. 2. Results of LOWTRAN7 simulations 
using five relative humidity profiles 
for TIMS Jasper Ridge data. 


Table 1. Surface Temperature 
Estimation using Split Window 
Technique 


[°C] 

Lake 

Track 

Lot 

TIMS observed 
average 

18.9 

40.7 

31.7 

(low) 

(17.6) 

(37.3) 

(29.6) 

(high) 

(19.6) 

(42.3) 

(33.1) 

Ch. 3 - Ch. 1 

1.6 

5.0 

2.0 

Ch. 5 - Ch. 6 

0.6 

1.8 

>0.1 

Ground Truth 

22.0 

48.0 

38.0 


Midlat. Summer 


average 

21.6 

51.7 

39.6 

(low) 

(21.2) 

(50.5) 

(37.3) 

(high) 

Radiosonde 

(22.0) 

(53.4) 

(43.3) 

average 

20.6 

48.2 

36.8 

(low) 

(20.3) 

(47.8) 

(35.6) 

(high) 

(20.9) 

(48.6) 

(39.1) 

Split Window 




Ch. I and 3 

21.0 

49.9 

34.1 

Ch. 5 and 6 

21.3 

47.9 

32.6 
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Figure 1 shows the estimated temperature of the average of 500 lines of TIMS data 
of the Pacific Ocean, from flight line 94, collected on September 30, 1988, at 1931 GMT. In 
Figure 1 A, with no atmospheric corrections, estimated temperature decreases away from 
nadir (the center of the scan line). In Figure IB, a LOWTRAN modeled correction (Kneizys, 
et al, 1985), using local radiosonde data and instrument scan angle information, results in 
reversed limb darkening effects for most bands, and does not adequately correct all bands to 
the same temperature. The atmosphere tends to re-radiate energy at the wavelengths at 
which it most absorbs, and thus the overall difference between corrected and uncorrected 
temperatures is approximately 4° C, despite the average LOWTRAN calculated 
transmittance of only 60% between 8. 1 and 1 1 .6 pm. An alternative approach to 
atmospheric correction is a black body normalization. This is done by calculating a 
normalization factor for each pixel position and wavelength, which when applied to Figure 
1 A, results in a single calculated temperature, as would be expected for a gray body with near 
uniform emittance. 

The black body adjustment is based on the atmospheric conditions over the sea. 
Figure 2 shows the ground elevation profile along the remaining 3520 scan lines 
(approximately 10 km) of flight line 94, up the slopes of Kilauea, determined from aircraft 
pressure and laser altimeter data. This flight line includes a large amount of vegetation that is 
clearly discernible on the radiance image, being much cooler than the surrounding rocks. For 
each of the 3520 scan lines, pixels were classified as vegetation or "other" (see Figure 3). A 
moving average of 51 lines was applied to the composite vegetation emittance (Warner and 
Levandowski, 1990) for each scan line, to reduce noise (Figure 4). Assuming vegetation to 
be like water, and to act as gray body with an emittance of 0.986 across the spectrum, Figure 
4A shows that the LOWTRAN induced artifacts are severe, and other than for the 0.9.9 pm 
channel, not significantly different from applying no corrections at all. As expected, with 
increasing elevation atmospheric effects are slightly reduced (Figure 4B), because moisture 
tends to be concentrated in the lowermost part of the atmosphere. The black body 
adjustment (Figure 4 C) is highly robust, and even at elevations nearly 600 meters above the 
sea, remains an alternative procedure for use in calculating emittance. 
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Figure 1 . Radiance data from the average of 500 lines of 1 IMS data of the Pacific Ocean, 
converted to estimated black body temperatures. A. Temperatures esitmated with no 
atmospheric corrections. B. Temperatures estimated using LOWTRAN modelled 
atmospheric parameters. 

Ground Elevation 



Figure 2. Nadir ground elevation for part of TIMS flightline 94, estimated from laser and 
pressure altimeter data. 


Number of pixels per scan line classified as vegetation 



Figure 3. Number of pixels in each scan line classified as vegetation. 
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Figure 4. Calculated six band emmitance for vegetation versus scan line number. Ground elevation increases with increasing line number (see Figure 2.) A. Lowtran corrected 
(Lowermost curve is 9.9 Jim). B. No corrections applied. (Lowermost curve is 8.4 Jim.) C. Black body adjustment 
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The effective flux incident upon the detectors of a thermal sensor, after it has been 
corrected for atmospheric effects, is a function of a non-linear combination of the 
emissivity of the target for that channel and the temperature of the target. The sensor 
system cannot separate the contribution from the emissivity and the temperature that 
constitute the flux value. 

In this paper, we describe a method that estimates the bounds on these temperatures 
and emissivities from thermal data. This method is then tested with remotely sensed data 
obtained from NASA’s Thermal Infrared Multispectral Scanner (TIMS) - a 6 channel 
thermal sensor. Since this is an under-determined set of equations i.e there are 7 
unknowns ( 6 emissivities and 1 temperature ) and 6 equations (corresponding to the 6 
channel fluxes), there exist theoretically an infinite combination of values of emissivities 
and temperature that can satisfy these equations. Using some realistic bounds on the 
emissivities , bounds on the temperature are calculated. These bounds on the temperature 
are refined to estimate a tighter bound on the emissivity of the source. An error analysis 
is also carried out to quantitatively determine the extent of uncertainty introduced in the 
estimate of these parameters. This method is useful only when a realistic set of bounds 
can be obtained for the emissivities of the data. In the case of water the lower and upper 
bounds were set at 0.97 and 1.00 respectively. 

Five flights were flown in succession at altitudes of 2 km (low), 6 km (mid), 12 km 
(high), and then back again at 6 km and 2 km. The area selected was the Ross Barnett 
reservoir near Jackson, Mississippi. The mission was flown during the predawn hours of 
Feb. 1, 1992. Radiosonde data was collected for that duration to profile the characteristics 
of the atmosphere. Ground truth temperatures using thermometers and radiometers were 
also obtained over an area of the reservoir. The results of two independent runs of the 
radiometer data averaged 7.03±.70 for the first run and 7.31±.88 for the second run. 

The results of the algorithm yield a temperature of 7.68 for the low altitude data to 
8.73 for the high altitude data. 
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Surface radiometric measurements 



Distribution of reservoir surface temperature 5AM 2/1/92 



Variation of Ground Temperature with altitude. 
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Average of brightness temperature in all channels. 


Flt.Ln. 

Mean 

Sigma 

Mean-sigma 

Mean+Sigma 

1 

6.51 



6.816 

2 

7.45 



7.728 

3 

7.86 


7.626 

8.094 

4 

8.38 

wmmmwm 

8.098 

8.662 

5 

7.39 

0.321 

7.069 

7.711 


Groung temperature using the bounds algorithm. 


Fit. Ln. 

Mean 

Sigma 

Mean-sigma 

Mean+Sigma 

1 

7.68 

0.422 

7.258 

8.102 

2 

8.74 


8.236 

9.244 

3 

8.73 


8.350 

9.110 

4 

9.29 

0.348 

8.942 

9.638 

5 

8.22 

0.410 

7.811 

8.630 
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Fractal geometry is increasingly becoming a useful tool for modeling natural 
phenomenon. As an alternative to Euclidean concepts, fractals allow for a more accurate 
representation of the nature of complexity in natural boundaries and surfaces. Since they 
are characterized by self-similarity, an ideal fractal surface is scale-independent; i.e at 
different scales a fractal surface looks the same. This is not exactly true for natural 
surfaces. When viewed at different spatial resolutions parts of natural surfaces look alike 
in a statistical manner and only for a limited range of scales. 

In this paper, images acquired by NASA’s Thermal Infrared Multispectral Scanner are 
used to compute the fractal dimension as a function of spatial resolution. Three methods 
are used to determine the fractal dimension - Shelberg’s line-divider method, the 
variogram method and the triangular prism method. A description of these methods and 
the result of applying these methods to a remotely-sensed image is also presented. 

Five flights were flown in succession at altitudes of 2 km (low), 6 km (mid), 12 km 
(high), and then back again at 6 km and 2 km. The area selected was the Ross Barnett 
reservoir near Jackson, Mississippi. The mission was flown during the predawn hours of 
Feb. 1, 1992. Radiosonde data was collected for that duration to profile the characteristics 
of the atmosphere. This corresponds to 3 different pixel sizes - 5m , 15m and 30m. 
After, simulating different spatial sampling intervals within the same image for each of 
the 3 image sets, the results are cross-correlated to compare the extent of detail and 
complexity that is obtained when data is taken at lower spatial intervals. 

introduction 

The advent of fractal analysis measures has been alluded to as one of the four most 
significant scientific concepts of the 20th century, with a scientific impact similar to that 
created by quantum mechanics, the general theory of relativity, and the development of the 
double-helix model in DNA structure. (Clarke et al. 1991) Since the development of the 
fractal concept by Mandelbrot (Mandelbrot 1967, Mandelbrot 1977), a number of methods 
for calculating fractal dimensions have been developed and applied to various spatial 
problems. (Goodchild 1980, Goodchild 1982, Burrough 1981, Mark et al. 1984, 
Goodchild et al. 1987, Krummel et al. 1987, Milne 1991, Lam et al. 1992) Fractals, 
however, have seen only limited employment for analysis of remote sensing data. 
(DcCola 1989, Lam 1990) As noted in Lam (1990), fractals offers significant potential 
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for improvement of measurement and analysis of spatially and spectrally complex remote 
sensing images. The fractal dimension of remote sensing data could yield quantitative 
insight on the spatial complexity and information content Thus, remote sensing data 
acquired from different sensors and at differing spatial and spectral resolutions could be 
compared and evaluated based on fractal measurements. The fractal dimensions derived 
from remote sensing data could also be compared with other measures of spatial 
complexity to better understand the significance of the spatial interrelationships present 
within image data. Moreover, if fractal dimensions are shown to be unique to different 
types of remote sensing data, these values could be used as control for the simulation of 
fractal surface generation of remotely sensed images. 

Outside of the potential offered by fractal analysis, the only other method for 
measuring the spatial frequency as an index of variability and complexity within a remote 
sensing image has been the two-dimensional Fourier transform technique (2D-FFT). In 
spite of being computationally intensive, the 2D-FFT technique has been successfully 
used for spatial processing of image data, such as filtering. (Moik 1980) The only 
information that a 2D-FFT can provide to classify data in a spatial manner is the 
spectrum of the data. For example, urban areas have a higher spatial frequency than rural 
areas. The information is very qualitative and does not readily lend itself to quantitative 
interpretation. With the emergence of fractals, it may be possible to lend that quantitative 
analysis of these spatial variations. 

This paper describes the adaptation and implementation of three methods that have 
been successfully applied to compute fractal dimensions from multiple scaled remote 
sensing data. These are the Shelberg or line-divider method, (Shelberg et al. 1982, 
Shelberg et al. 1983) the triangular prism method, (Clarke et al. 1991), and the variogram 
method. (Mark et al. 1984) These techniques have been implemented on a self-contained 
menu-driven PC-compatible image interpretation software package (Note: This package is 
also available for the UNIX workstation environment, without the menu structure, as a 
set of routines). This interactive program, written in 'C', allows the user to analyze their 
results without having to spend considerable effort in programming these methods for 
fractal computation. The results are accessible to the user on the screen as well as from 
an ASCII File for future use in graphing and more intensive interpretation. Also, this PC- 
compatible software utilizes almost no special purpose hardware (except for a VGA 
monitor), which permits wide distribution to other users. 

On entering the program, the user encounters an image on the screen with a menu 
adjacent to it. The menu allows the user to perform analyses using either of the three 
fractal computation methods. The fractal dimension is computed for a user- specified 
region within the image. This region is outlined as a rectangular box on the image. The 
user can then interactively move the box anywhere or change the size of the box. Also, 
some of these methods allow for analysis to be performed using different values for 
specific internal parameters. These parameters can also be changed interactively or can be 
preset by the user. The following is a discussion of the three fractal computation 
methods, their implementation using the program, and a description of the results 
obtained from each fractal calculation method. 
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1, G EOLO G IC SE TTING 

The Iron Hill carbonatite-alkalic igneous rock complex is in the Powderhom mining 
district, approximately 40 km south -south west of Gunnison, Colorado . The complex, which 
occupies about 30 km 2 , was emplaced in metasedi memary and metavolcanic rocks during the 
late Precambrian or early Cambrian (Olson and Hedlund, 1981). The main rock types in the 
complex, from oldest to youngest, are fenite, pyroxenite, uncompahgrite, ijolite, nepheline 
syenite, and doiomidc carbonarite. The carbonadte is Iimonidc and forms an ellipdcally 
shaped 4km 2 stock. Calcidc and doiomidc carbonadte dikes are also numerous throughout the 
complex and in the pre-existing rocks. Pyroxenite is the most widespread rock type within 
the complex, but pyroxene is extensively altered to biodte, phlogopite, and vermiculite. 

Fenite, which formed through Na, K-metasomadsm of the country rocks (Heinrich, 1980), 
typically contains more feldspar and less quartz than the equivalent unaltered country rocks. 
The other alkali c rock types axe less widespread and less well exposed. Parts of the complex 
are covered by Oligocene ash-flow tuff and alluvial, colluvial, and glacial deposits. 

Sagebrush and grass cover is moderately dense to very dense at low to intermediate 
elevations; coniferous tree cover is dense at high elevadons and on some north-facing slopes 
at lower elevadons. 

2. DATA PROCESSING 

TIMS data were acquired of the Iron Hill area on July 22, 1990 from a mean aidtude 
above terrain of 4000m. Weather conditions prior to the flight indicated that the ground was 
dry. Initial data calibration was conducted by assuming a linear relationship berween the 
observed blackbody signals and the band pass fluxes at the hot and cold blackbody 
temperatures convolved with the spectral responses for the six channels (PaLluconi and 
Meeks, 1985). The coefficients of this linear relationship were then used to convert the 
observed signals to equivalent fluxes. A generic spectral atmospheric correction for the path 
radiance, transmission, and downward sky radiance based on MODTRAN (Berk et al., 1989) 
was then employed. 

Several methods have been used for extracting spectral emissivity information from 
TIMS data (Gillespie et al., 1986; Kahle and Rowan, 1980; Lahren et aJL, 1988; Watson and 
Raines, 1989; Watson et al., 1990; Gillespie, 1986; Hook et al., 1990; Kealy and 

Gabell, 1990). Decorrelation stretches have been commonly used (Gillespie et al., 1986), 
although they are of limited value because they do not yield interpretable spectral 
information. 

In this study a new algorithm (Watson, 1992a) is used to compute spectral emissivity 
ratios, independent of any emissiviry assumptions. This algorithm has the advantage that any 
of the possible emissivity ratios can be computed and, thus, a large variety of composite ratio 
images can be constructed, which permits examination of various geologic hypotheses based 
on the spectral properties of the surface materials. 

Prior to processing the TIMS data, evaluation of laboratory reflectance spectra of 
field samples identified 13 TIMS ratios that might be useful for distinguishing the main rock 
types. 
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Initial analysis of these ratio images resulted in selection of the following images for 
detailed evaluation: 2 J U J /,, A / u *l u V 6 , 3 / 6 , and 2 / 5 . A vegetation mask, which was based on 
high 3 / 2 digital numbers, was used to reduce the influence of vegetation on the contrast 
stretches and, hence, increase lithologic discrimination. The cause of the high digital 
numbers in vegetated areas in this and the other ratio images is not clear. The ratio images 
were then spatially filtered to remove scanline noise using a two-dimensional fast-Fourier 
transformation (Watson, 1992b) and registered to a 1:24,000 topographic map. The 
registration was accomplished by first using a small set of control points to perform an affine 
transformation and, then, an inverse distance weighed algorithm (Hummer-Miller, 1990) was 
employed using over 100 control points. 

3, INTEKPRETATIQN and conclusions 

Most of the alkalic igneous rocks (pyroxenite, ijolite, and uncompahgrite) are readily 
distinguished from the silica-saturated rocks (quartzite, granite, and ash-flow tuff)in each of 
several color-ratio composite images, especially those including the 2 / w J / lf l / 4 , 3 / 4 , and 3 / 1 
ratios. This implied large spectral emissivity contrast is consistent with the markedly 
different laboratory reflectance spectra of these two rock groups. In general, quartzite, 
granite, and ash-flow tuff are spatially separable in the color-ratio composite images, but 
color variations within the alkalic-rock group may reflect the proportions of rock, soil, and 
vegetation, as well as rock composition differences. Although the carbonatite stock is 
distinctive from the alkalic rocks, it is difficult to distinguish from some of the other rock 
types. 


Distinction of the carbonatite from tuff and from some granitic areas, and separation 
of the fenite from granite were only feasible using color-ratio composite images which had 
had the vege&tion mask applied prior to contrast stretching. Several different masked color- 
ratio composite images were needed to display subtle spectral emissivity difference that are 
critical for mapping these lithologic variations. Discrimination of the carbonatite is 
hampered by the location of the weak CO, reststrahlen feature between TIMS channels 5 and 
6 . 
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INTRODUCTION 

S-type and p-type pahoehoe record different mechanisms and vigors of activity 
within an active flow field (e.g. Swanson 1973; Hon 1991; Multhaup & Walker 1992). 
There is some controversy about what these mechanisms are exactly, and this study was 
undertaken with the idea that an accurate map of the two surfaces within a pahoehoe 
flow field could be helpful in solving the problem. TIMS allows discrimination 
between s-type and p-type pahoehoe, and this ability has been used to map the two 
surface types on the Mauna Iki satellitic shield (southwest rift zone, Kilauea Volcano). 

TIMS has previously been used to discriminate a'a from pahoehoe as well as to 
determine relative age relationships of different flows (e.g. Kahle et al. 1988; Abrams 
et al. 1991). Although inter-flow variation was minor in the data published by these 
authors, a second goal of the work presented here is to understand such variations to 
better constrain intra-flow differences used for age dating. 

DATA USED 

The TIMS data were collected in November 1985 and have a spatial resolution of 
* 10 m. An image consisting of the first 3 principal components (PC) and detailed field 
checking were used to derive a unit map (Figure 1). The PC image was useful for unit 
discrimination but because we found that the spectral character of the two pahoehoe 
surface types changed downflow, no rigorous classification scheme was devised. 

After constructing the unit map, we extracted broad-band TIMS spectra to aid in 
explaining the mechanism behind the ability to discriminate the two lava types. We 
also made measurements of surface temperature in the field using a thermocouple. 
Quantitative analysis of the PC image, the spectra, and the surface temperature 
measurements (in collaboration with colleagues at JPL) is ongoing. 

S-TYPE AND P-TYPE PAHOEHOE 

S-type and p-type pahoehoe differ both in their character and their distributions 
within a flow field. In short, s-type pahoehoe is vesicle-rich (s stands for spongy), and 
is erupted directly from a lava tube. The surface of s-type consists of stretched and 
broken vesicles, and often spalls off shortly after emplacement. P-type pahoehoe is 
vesicle poor, and its chilled margin resembles obsidian. The reason for the paucity of 
vesicles comprises the present controversy. The most commonly-held icjea is that p- 
type pahoehoe resides within a flow field for a period of time long enough to degas. If 
that particular storage area is then invaded by fresher lava, the degassed lava is pushed 
out on to the surface. This explains the transition from early p-type to later s-type lava 
often seen while watching a lava breakout. A large amount of degassing is known to 
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take place during flow in tubes (also recently investigated using TIMS: Realmuto et aL 
1991). 

The competing idea is that p-type pahoehoe is actually richer in gas than s-type, 
and the pressure provided by the overlying flow field in which it resides causes solution 
of gas back into the lava. P-type pahoehoe often has a lower aspect ratio, possibly 
reflecting a lower viscosity generated by the higher dissolved gas content. Arguing 
against this is the fact that reabsorbtion of gas is endothermic whereas p-type pahoehoe 
is often hotter than s-type upon eruption. 

S-TYPE AND P-TYPE DISTRIBUTION AT MAUNA IKI 

In Figure 1, the "main shield" and "downrift flow" were determined by field 
mapping and studying the eruption chronology. Within these two zones, the 
distribution and relative amounts of s-type and p-type pahoehoe differ greatly. On the 
main shield, p-type pahoehoe is restricted to 3 small patches along the northwest 
margin of the flow. It is here that lava flows from the Mauna Iki summit ponded 
against pre-existing topography. This is a condition conducive to stagnation of the flow 
and consequent inflation. The p-type pahoehoe erupted out of storage within this 
stagnated flow. The major activity on the main shield was the overflowing of relatively 
gas-rich lava from a small lava pond, and the formation of p-type pahoehoe was an 
exception. 

The downrift flow was emplaced by lava tubes, and the distribution and relative 
abundance of the two pahoehoe types are very different. They make up roughly equal 
proportions of the flow and at contacts p-type overlies s-type. Overall, the Mauna Iki 
eruption had a low eruption rate, and this was particularly true of the final lava tube 
stage. This low eruption rate meant that once the flow field was established, it was 
difficult to supply all of it with fresh lava. 

The distribution of p-type pahoehoe on the flow field can be used to constrain the 
conditions required for formation. First, there is a paucity of p-type pahoehoe where 
lava tubes were well-established (i.e. on the shield and the main axis of the downrift 
flow). This indicates that initially, lava erupted directly from tubes was too gas-rich to 
form p-type. The downrift flow field averages "1 m in thickness, and is even thinner 
along its margins. The large amount of p-type pahoehoe compared to the flow volume 
as a whole argues against the existence of an overburden capable of pressurizing the 
lava. The field relations suggest instead, that even though gas-rich lava established the 
distal end of the flow initially, in the waning later stages all the lava degassed during 
sluggish travel down the tubes prior to being erupted onto the surface. 

Kipukas (inliers) of s-type tumuli form semi-linear trains branching off from the 
main axis of the downrift flow. These may indicate the lines of subsidiary early-formed 
tubes, and the ability to map tubes long after a lava flow has stopped erupting would be 
a very useful tool. We are in the process of conducting this same type of analysis on 
the Mauna Ulu flow field where the positions of lava tubes were closely monitored 
during the eruption. 

FUTURE WORK AND CONCLUSIONS 


The TIMS data indicate that at Mauna Iki the radiometric temperature of s-type 
pahoehoe is hotter than that of p-type, supported by the fact that PCI discriminates the 
two types most strongly. Further analysis of the TIMS data, as well as more complete 
field temperature measurements, will help to determine if these differences are due to 
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different physical temperatures or emissivities. In conclusion, it is possible to 
discriminate pahoehoe surface types using TIMS data. Quantitative analyses of the 
causative relationships, however, have not been completed, and these are the next 
goals. With regard to the formation of p-type pahoehoe, the large area covered by p- 
type, including areas where the flow as a whole is thin, suggests that the availability of 
an overburden is not required for its formation, in turn supporting the degassed origin. 

I am very grateful to Harold Garbeil for helpful discussions and for his image 
processing skills. 
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Figure 1. Unit map derived from TIMS data for the Mauna Iki flow. 
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L INTRODUCTION 


Accurate interpretation of thermal infrared data depends upon the understanding 
and removal of complicating effects. These effects may include physical mixing of var- 
ious mineralogies and particle sizes, atmospheric absorption and emission, surficial coat- 
ings, geometry effects and differential surface temperatures. The focus of this study is the 
examination of the linear spectral mixing of individual mineral or endmember spectra. 
Linear addition of spectra, for particles larger than the wavelength (SaUsbury et al, 1987), 
allows for a straight-forward method of deconvolving the observed spectra, predicting a 
volume percent of each endmember. The ’forward analysis' of linear mixing (comparing 
the spectra of physical mixtures to numerical mixtures) has received much attention 
(Thomson & Salisbury 1991; Christensen et al, 1986). The reverse approach of un-mixing 
thermal emission spectra has been examined with remotely sensed data (Adams et al, 

1989; Gillespie et al, 1990), but no laboratory verification exists. Understanding of the 
effects of spectral mixing on high resolution laboratory spectra allows for the extrap- 
olation to lower resolution, and often more complicated, remotely gathered data. 

'niermal Infrared Multispectral Scanner (TIMS) data for Meteor Crater, Arizona 
were acquiredTn September, 1987. The spectral un-mixing of these data gives a unique 
test of the laboratory results. Meteor Crater (1.2 km in diameter and 180 m deep) is 
located in north-central Arizona, west of Canyon Diablo (Shoemaker & Kieffer, 1974). 

The arid environment, paucity of vegetation and low relief make the region ideal for 
remote data acquisition. Within the horizontal sedimentary sequence that forms the upper 
Colorado Plateau, the oldest unit sampled by the impact crater was the Permian Coconino 
Sandstone. A thin bed of the Toroweap Formation, also of Permian age, conformably 
overlays the Coconino. Above the Toroweap lies the Permian Kiabab Limestone which, 
in turn, is covered by a thin veneer of the Moenkopi Formation, The Moenkopi is 
Triassic in age and has two distinct sub-units in the vicinity of the crater. The lower, 
Wupatki member, is a fine-grained sandstone, while the upper, Moqui member, is a 
fissile siltston z (Shoemaker & Kieffer, 1974). Eje cta fro m these units are preserved as 
inverted stratigraphy up to 2 crater radii from the rim. The mineralogical contrast LL. 
between the units, relative lack of post-empiacement erosion (Grant & Schultz, 1989) and 
ejecta mixing provide a unique site to apply the un-mixing model. Selection of the 
aforementioned units as endmembers reveals distinct patterns in the ejecta of the crater. 


2 . Laboratory Spectra 


34 


li,l !i liUiiliiiiiliiiiliUilJJii 



Particles between 355 and 500|im were used to reduce the effects of volume scat- 
tering; the choice of endmember minerals approximates the mixing of Kiabab Limestone 
and Coconino Sandstone in the crater ejecta. Laboratory thermal emission spectra (Fig. 1) 
of the powders were acquired on a Mattson Cygnus 100 FTIR interferometer/spectro- 
meter. Absolute emissivity was obtained using the technique described by (Christensen & 
Harrison, 1992). The endmember components, mixed by volume into binary mixtures, 
were agitated after each of the five spectral acquisitions. This process assured accurate 
spectral sampling and charted the daily variation of the spectrometer, which varied by less 
than 1% absolute emissivity. The five runs were then averaged to produce a final 
spectrum. The calcite endmember spectra shows the characteristic absorption band 
associated with the vibration of the carbonate ion (® 1550 cm 1 ) and the large reststrahlen 
band produced by the Si-0 bending mode of quartz (-1150 cm" 1 ) (Salisbury et al , 1987). 
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FIG. 1. Emission 
spectra of 355-500|j.m 
quartz and calcite. The 
dashed line shows the 
physical mixture of 40% 
quartz, 60% calcite. 
Errors (Ac) in the model 
fit for the 40% quartz 
mixture are shown in the 
bottom curves. The solid 
line is the two endmem- 
ber case, while the dotted 
line indicates the 
improved error after 
blackbody correction. 


The spectra of the two endmembers were used as reference inputs to the model. A 
least-squares fit of the data results in the percentage of each endmember as well as the 
"goodness-of-fif of the model to the data (rms error) (Adams et al, 1989; Gillespie et al, 

1990). Figure 2 shows the results for the quartz endmember. Predicted percentages were 
within 8 % of actual values and validate the linear assumption for particles larger than the 
wavelength. The largest rms errors (Fig. 1) occur over the strong absorption features, and 
are due to a band shallowing in the spectra of the physical mixtures. The behavior over 
these regions, a function of the photon path length and the imaginary part of the index of 
refraction k(X), appears to be slightly non-linear. The use of unimodal grain sizes virtually 
eliminates a change in the path length, indicating the non-linear behavior must be related 
to the change in the overall k(X) of the mixture. The band shallowing can be approximated 
by the addition of blackbody emissivity (Ramsey &. Christensen , 1992). 



FIG. 2. Model 
predicted results for the 
percentage of quartz. 
Squares indicate the 
values obtained using 
two endmembers; the 
diamonds show the 
improved results using 
blackbody emissivity as 
a third endmember. 
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Using a constant emissivity of 1 .0 as a third endmember improved the results to better 
than 5%, reducing the errors over the bands. The model, re-run on the lab spectra 
convolved to TIMS resolution, predicted the percentages to within 10 %; however, no 
attempt was made to model the effect of additional noise. 


3. Meteor Crater Image analysis 


Standard image enhancement techniques such as the decorrelation stretch do not 
account for large temperature differentials within the scene and tend to stretch microphonic 
noise as well. A spectral un-mixing analysis eliminates these problems, however. For 
this study, six emissivity images were extracted from the TIMS data using the spectrum 
normalization technique (Realmuto, 1990). An improvement in the algorithm allows for 
interactive linear stretching of the emissivity images for a maximum spectral contrast prior 
to input into the un-mixing program. The aforementioned stratigraphic units were used as 
endmembers for the model (5 total). These endmembers accurately fit the data to within ± 
2 DN except for areas over the buildings and road. The Coconino endmember image 
clearly reveals a NE trending windstreak (Grant &. Schultz , 1989), which is a thin (<50 cm) 
veneer of eolian-derived material. In addition, subtle patterns within the Kiabab indicate 
the true extent of the ejecta. This supports the interpretation of (Grant & Schultz, 1989) 
who state that the ejecta is well-preserved below the thin alluvial mantling deposits derived 
from Holocene erosion. These patterns are not evident in the decorrelation stretched image 
and must be field verified. Differentiation of slight mineralogical changes within the 
Moenkopi are also clearly distinguished and allow for improved eject mapping. 
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TIMS data were acquired over Mount Etna, Sicily, on July 29, 1986 (Bianchi et 
al. 1990). The volcanic activity at that time was characterized by the steady effusion of 
gas from the Bocca Nuova (BN), Chasm, and Southeast craters. The Northeast crater, 
quiet at the time of the TIMS overflight, was the site of Slrombolian eruptive activity 
between July 31 and September 24, 1986 (SEAN 1986). 

In aerial photographs of the Etna summit region acquired during the TIMS 
overflight, the S 02 -rich plume is visible due to the scattering of sunlight by the entrained 
aerosols. In the TIMS imagery, the plume is revealed by the strong absorption of S0 2 
between 8 and 9 pm. This absorption feature falls within the first three channels of TIMS, 
with the strongest absorption falling within Channel 2. Following decorrelation 
processing, the plume is visible in color-composites of TIMS channels 2, 3, and 5. 

To estimate the concentration of S0 2 within the plume, the LOWTRAN 7 
radiative transfer code was used to model to radiance spectra measured by TIMS. Inputs 
to LOWTRAN included ground temperature and cmissivity, vertical profiles of 
atmospheric temperature and humidity, and the vertical position of S0 2 within the plume. 
The ground temperatures and cmissivity spectra were derived from the TIMS 
measurements through the use of a curve-fitting algorithm (Warner and Levandowski 
1990, Realmuto 1990). Emissivity spectra recovered outside of the plume were applied 
to the ground beneath the plume where appropriate. The ground temperature estimates 
were modified interactively by fitting LOWTRAN-gcncratcd radiance values to those 
measured in TIMS channels 4, 5, and 6, the channels that were not strongly affected by 
the presence of S0 2 in the atmosphere. The atmospheric pressure, temperature, and 
humidity profiles were recorded by radiosondes launched from Trapani, which is located 
approximately 220 km west of Mount Etna. 

To evaluate the accuracy and precision of the estimation procedure, S0 2 
concentrations were recovered over the same piece of ground using the data from three 
successive TIMS overflights. The flights were separated by intervals of approximately 15 
minutes. Assuming that the gas output of the volcano was constant during the course of 
the flights, the concentration of S0 2 in the chosen portion of the plume was 26.3 ± 4.1 
mg/m . This estimate compares favorably with S0 2 concentrations of 22.84 ± 2.3 and 
57.1 ± 5.7 mg/m 3 measured near the rim of BN crater in 1978 and 1979, respectively 
(Jaeschke et al. 1982). Assuming a plume width of 1 km and a wind speed of 5 m/s, the 
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S0 2 concentrations derived from the TIMS data correspond to flux rales ranging between 
5400 and 7600 tons/ day. These S0 2 flux rates arc within the range of published flux 
rates from Mount Etna (cf. Haulet et al. 1977, Zcttwoog and Haulet 1978, Malinconico 
1979), but the TIMS-dcrivcd (lux rates arc higher than the average flux rates published in 
the literature. 


The sensitivity of the estimation procedure to variations in the inputs to the 
LOWTRAN code was also evaluated. The estimation procedure was most sensitive to 
ground temperature; a 10 degree (Celsius) variation in ground temperature translated into 
a 50% change in the estimation of S0 2 concentration. Changes in the vertical distribution 
of S0 2 within the plume resulted in 20% changes in the estimates of concentration. 
Replacing the ground emissivity with a uniform value of 0.95 resulted in a 10% change 
in the estimates of concentration. The substitution of a climatological atmospheric profile 
(mid-latitude summer) for thcTrapani radiosonde measurements resulted in a 20% change 
in the estimates of S0 2 concentration. 

The TIMS data provided a near- instantaneous measurement of S0 2 
concentration of the plume in the summit region of Mount Etna. Data such as these could 
aid in the study of temporal and spatial variations in the S0 2 content of the plume. It is 
recommended that future TIMS flights over venting volcanoes be accompanied by 
ground support crews. These personnel could collect ground temperature and emissivity 
measurements as well as launch radiosonde balloons. Despite the fact that the Mount 
Etna TIMS mission was not supported by a ground crew, the S0 2 concentrations derived 
from the TIMS data were in good agreement with published values. The availability of 
ancillary data can only improve the accuracy of the TIMS estimation procedure. 
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The Advanced Sensor Development Laboratory (ASDL) at the Stennis Space Center 
develops, maintains and calibrates remote sensing instruments for the National 
Aeronautics & Space Administration. To perform system design trade-offs, analysis, and 
establish system parameters, ASDL has developed a software package for analytical 
simulation of sensor systems. This package has been useful in modeling and analyzing 
the performance of various sensors such as the TIMS. This package called "Analytical 
Tools for Thermal InfraRed Engineering " - ATTIRE simulates the various components 
of a sensor system. The software allows each subsystem of the sensor to be analyzed 
independently for its performance. These performance parameters are then integrated to 
obtain system level information such as SNR, NER, NETD etc. This paper describes the 
uses of the package and the physics that were used to derive the performance parameters. 

In addition, ATTIRE can be used as a tutorial for understanding the distribution of 
thermal flux or solar irradiance over selected bandwidihs of the spectrum. This spectrally 
distributed incident flux can then be analyzed as it propagates through the subsystems that 
constitute the entire sensor. ATTIRE provides a variety of functions ranging from 
plotting black-body curves for varying bandwidihs and computing the integral flux, to 
performing transfer function analysis of the sensor system. 

The package runs from a menu-driven interface in a PC-DOS environment. Each 
sub-system of the sensor is represented by windows and icons. A user-friendly mouse- 
controlled point-and-click interface allows the user to simulate various aspects of a sensor. 

The package can simulate a theoretical sensor system. Trade-off studies can be easily 
done by changing the appropriate parameters and monitoring the effect on the system 
performance. The package can provide plots of system performance versus any system 
parameter. A parameter (such as the entrance aperture of the optics) could be varied and 
its effect on another parameter (e.g., NETD) can be plotted. A third parameter (e.g., the 
obscuration) could be varied for each plot and several plots obtained on the same graph. 
The menu for such " Y vs X plots for different values of Z ” currently contains various 
such options. The package also allows the user to create customized work-sheets of the 
simulated system and save the analysis for interface with or retrieval to other packages. 

The emissivity, atmospheric transmission and the optical transmission default as 
constants over the specified spectral bandwidihs. There is an option for making these 
parameters spectrally variable. If more than one of the above-mentioned three parameters 
are spectrally variable, then it is possible that the upper and lower wavelength values as 
well as the resolution of the wavelength array may not be the same for all three arrays. 
The package performs an interpolation of the data to smooth out the curves and then 
projects them onto a common wavelength array for all the parameters. 
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Overview 


ATTIRE provides functions ranging from plotting blackbody curves and computing 
the integral flux, to performing radiation transfer function analysis and calculation of 
system performance parameters. 

The package runs from a menu-driven interface in a PC -DOS environment. On 
executing the software, the user encounters a screen with a resizeable window and icons. 
The window contains the channels of the current multispectral sensor being analyzed. 
Each icon represents a sensor subsystem such as the optics, source flux, atmospherics etc. 
ATTIRE has a mouse-driven control for resizing or moving the windows within the 
screen. The mouse control also allows for the point-and-click feature making the package 
more user-friendly. 

A horizontal menu within the window allows the user to change the parameters of 
the system. This menu contains items that allow the user to change the channel 
configuration, perform file I/O and do plotting operations. The number of channels is 
programmable through the horizontal menu item called CHAN. The bandwidths of each 
of these channels can also be varied from this menu-item. The entire sensor system can 
be saved as a worksheet in a file. This operation is done from the horizontal menu-item 
called FILE. The user can either retrieve an existing sensor system, save the current 
system or just browse through a system saved on disk before deciding to retrieve it 

The source emissivity, atmospheric transmission, and optical transmission are 
initially assumed constant over the specified spectral bandwidths. There is an option for 
making these parameters spectrally variable. If more than one of the above-mentioned 
three parameters are spectrally variable, then it is possible that the upper and lower 
wavelength values and the resolution of the wavelength array may not be the same for all 
three arrays. The software package performs an interpolation of the data to smooth out 
the curves and then projects them onto a common wavelength array. 

The spectrally distributed incident flux can be analyzed as it propagates through the 
subsystems. Fig. 1 shows the signal flow and the various parameters that affect the 
performance of the system. Each of the icons on the main menu represent one of these 
sub-systems. 

Summary 

The design of a visible through thermal IR sensor system requires a detailed analysis 
of how the input signal propagates through the system. The major components in 
developing a model for the sensor system are the source, atmosphere, optics, 
detector, spatial parameters, and preamplifier electronics. The final goal of 
the analysis is to determine the NER for the various spectral channels of a sensor system. 

In this paper, a simulation package "ATTIRE” for analyzing sensor systems was 
introduced. The package runs in a PC-DOS environment and consists of one executable 
program and several supporting files. The entire package fits on one high-density floppy 
disk. 


ATTIRE is a useful tool for performing design tradeoffs as it interrelates several 
aspects of the sensor system to yield performance parameters. It is also useful as a tool 
I for the understanding of the concepts of radiometry. 

i 
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Fig. 1 


Wavelength vs E Flu* 
microns v» Wem*-2 



Fig. 2 Solar Curve used in ATTIRE 



Fig. 3 Output from ATTIRE 
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1. DATA COMPILATION 

In the last several years, the NASA aircraft program has amassed 
a considerable set of remote sensing data. With the upcoming launch of 
the EOS instrument platform at the end of this decade, some type of mass 
data organization is necessary in order to efficiently archive both pre- and 
syn-EOS data. With this in mind, we have embarked on a series of 
Compact Disc volumes that will consist of data compilations for various 
areas. The first set in this series will be for Kilauea volcano on the island 
of Hawaii. Kilauea volcano makes up the southeast portion of the big 
island, and has been erupting along the east rift zone almost continuously 
since 1983. 

Since 1985, the NASA C- 1 30 aircraft has been deployed four times 
to the island of Hawaii. For each of these deployments, the complement 
of instruments included Zeiss metric cameras and the NS001 (Thematic 
Mapper Simulator) scanner, and all but the 1987 deployment also included 
the Thermal Infrared Multispectral Scanner (TIMS). Due to space 
restrictions, a single flight line of photographs for each area on Kilauea 
will be included on the disc, except in areas that underwent visible changes 
between deployments (i.e., the Pu’u O’o lava flow field). In addition to 
these C-130 data, the NASA DC 8 carrying the NASA/JPL Airborne 
multipolarization, multifrequency Synthetic Aperture Radar (AIRSAR) was 
deployed to Hawaii in 1990, mapping a large portion of Kilauea. The 
Compact Disc compilation will contain most of the useable data collected 
by these instruments over much of the volcano including multiple coverage 
of the active areas. 

To aid in the use of the aircraft data, several other ancillary data 
types will also be included on the disc. These include a Digital Elevation 
Model (DEM) of Kilauea, and several maps showing the progression of 
lava flows over time provided courtesy of the U.S. Geological 
Survey/Hawaii Volcano Observatory. We will also include ‘browse’ 
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images of a Thematic Mapper scene and one or two NOAA AVHRR 
passes of the big island. Infrared laboratory spectra for many samples 
taken from various locations on Kilauea will also be included. The 
samples have been collected over several years during many different field 
campaigns. 

2. DATA ORGANIZATION 

The CDRom volume set will contain public domain image display 
software that is compatible with Sun Unix, PC and Macintosh operating 
systems. ReadMe files will be included for each data type explaining the 
data storage format as well as describing its possible uses. For some of 
the data types (i.e., TIMS), common data manipulation processes will be 
explained, and in a few cases, the algorithms will be provided. 

The first disc in the set will contain examples of each of the data 
types. It will be designated as a sampler disc containing thoroughly 
processed versions of all data types covering two selected areas on the 
volcano. For these areas full color air photos and contrast enhanced NS001, 
decorrelation stretched TIMS, and fully processed SAR images will be 
provided for easy viewing. 

The CDRom package is intended to provide a consolidated package 
of a sizable portion of the NASA aircraft acquired data for Kilauea 
volcano. The aircraft data, as well as the other types of ancillary data, will 
all be accessible in an easy to use format, in a single place. It is hoped 
that researchers from many disciplines working on Kilauea volcano will 
find this compilation useful. 
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